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Abstract: We extend the recent work on fluid-gravity correspondence to charged black-branes by de- 
termining the metric duals to arbitrary charged fluid configuration up to second order in the boundary 
derivative expansion. We also derive the energy-momentum tensor and the charge current for these 
configurations up to second order in the boundary derivative expansion. We find a new term in the 
charge current when there is a bulk Chern-Simons interaction thus resolving an earlier discrepancy 
between thermodynamics of charged rotating black holes and boundary hydrodynamics. We have also 
confirmed that all our expressions are covariant under boundary Weyl-transformations as expected. 
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1. Introduction 

Modern theoretical physics provides mathematically precise descriptions of a bewildering variety of 
phenomena. The diversity of the phenomena studies by physicists encourages specialization and a 
consequent divergence in the field. Wc find it satisfying, however, that this push towards divergence is 
partially counterbalanced by periodic theoretical discoveries that unify - i.e. discover precise mathe- 
matical connections between - distinct fields of physics. In line with this tradition, recent string theory 
inspired studies of classical gravitational dynamics have found a precise mathematical connection be- 
tween a long distance limit of the Einstein equations of gravity and the Navicr Stokes equations of 
fluid dynamics. More specifically, it has recently been demonstrated that a class of long distance, 
regular, locally asymptotically AdSd+i solutions to Einstein's equations with a negative cosmological 
constant is in one to one correspondence with solutions to the charge free Navier Stokes equations in 
d dimensions [1-6,8]. ^ 

^There exists a large literature in deriving linearise hydrodynamics from AdS/CFT. See( [9] - [41]). There have been 
some recent work on hydrodynamics with higher derivative corrections [43,44]. 
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The connection between the equations of gravity and fluid dynamics, described above, was demon- 
strated essentiafly by use of the method of collective coordinates. The authors of [1,3-5,8] noted that 

there exists a d parameter set of exact, asymptotically AdSd+i black brane solutions of the gravity 
equations parameterized by temperature and velocity. They then used the 'Goldstone' philosophy to 
promote temperatures and velocities to fields. The Navier Stokes equations turn out to be the effective 
'chiral Lagrangian equations' of the temperature and velocity collective fields. 

This initially surprising connection between gravity in + 1 dimensions and fluid dynamics in d 
dimensions is beautifully explained by the AdS/CFT correspondence. Recall that a particular large N 
and strong coupling limit of that correspondence relates the dynamics of a classical gravitational the- 
ory (a two derivative theory of gravity interacting with other fields) on AdS^+i space to the dynamics 
of a strongly coupled conformal field theory in d flat dimensions. Now the dynamics of a conformal 
field theory, at length scales long compared to an effective mean free path (more accurately an equili- 
bration length scale) is expected to be well described by the Navier Stokes equations. Consequently, 
the connection between long wavelength solutions of gravity and the equations of fluid dynamics - 
directly derived in [1] - is a natural prediction of the AdS/CFT correspondence. Using the AdS/CFT 
correspondence, the stress tensor as a function of velocities and temperatures obtained above from 
gravity may be interpreted as the stress tensor fluid of the dual boundary field theory in its deconfined 
phase. 

Now consider a conformal field theory that has a conserved charge Q in addition to energy and 
momentum. This is especially an interesting extension of the hydrodynamics of the uncharged fluids 
since the hydrodynamics of many real fluids has a global conserved charge which is often just the 
number of particles that make up the fluid. The long distance dynamics of such a system is expected 
to be determined by the augmented Navier Stokes equations; V^T^^" = together with V ^Jq = 0, 
where the stress tensor and charge current are now given as functions of the temperature, velocity 
and charge density, expanded to a given order in the derivative expansion. The bulk dual description 
of a field theory with a conserved charge always includes a Maxwell field propagating in the bulk. 
Consequently the AdS/CFT correspondence suggests asymptotically AdS long wavelength solutions 
of appropriate modifications of the the Einstein Maxwell equation are in one to one correspondence 
with solutions of the augmented Navier Stokes equations described above. 

This expectation of the previous paragraph also fits well with the collective coordinate intuition 
described above. Recall that the Einstein Maxwell equations have a well known d-\-l dimensional set 
of charged black brane solutions, parameterized by the brane temperature, charge density and velocity. 
It seems plausible that the effective Goldstone equations, that arise from the promotion of these d+1 
dimensional parameters to fields, are simply the augmented Navier Stokes equations. In this paper 
we verify the expectations via a direct analysis of the relevant bulk equations. More concretely, we 
generalize the work out in [1] to set up a perturbative scheme to generate long wavelength solutions 
of the Einstein Maxwell equations plus a Chern Simons term (see below for more details) order by 
order in the derivative expansion. We also implement this expansion to second order, and thereby find 
explicit expressions for the stress tensor and charge current of our dual fiuid to second order in the 
derivative expansion. 

In this paper we work with the Einstein Maxwell equations augmented by a Chern Simon's term. 
This is because the equations of IIB SUGRA on AdSs x S^ (which is conjectured to be dual to A/" = 4 
Yang Mills) with the restriction of equal charges for the three natural Cartans, admit a consistent 
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truncation to this system. Under this truncation, wc get the following action : 



S = 



1 




(1.1) 



The value of the parameter «; for = 4 Yang Mills is given by k = —1/(2^3) - however, with a view 
to other potential applications we leave k as a free parameter in all the calculations below. Note in 
particular that our bulk Lagrangian reduces to the true Einstein Maxwell system at k = 0. 

Our expressions for the charge current and the stress tensor of the fluid are complicated, and 
are listed in detail in (). We would however to point out an important qualitative feature of our 
result. Already at first order, and at nonzero k, the charge current includes a term proportional 
to = e'"'^"u^V^?iA. The presence of this term in the current resolves an apparent mismatch 
between the predictions of fluid dynamics and the explicit form of charged rotating black holes in 
IIB supergravity reported in [42] . Note that despite the presence of the e symbol, this term is parity 
symmetric, when accompanied by a flip in the charge of the brane, and so is CP symmetric in agreement 
with expectations of the CP symmetry oi Af — 4 Yang Mills theory. 

As we have explained above, the reduction of boundary field theory dynamics is expected to 
reduce to field theory dynamics only at long wavelength compared to an effective mean free path or 
equilibration length scale. All the; gravitational constructions of this paper also work only in the same 
limit. It is consequently of interest to know the functional form of the equilibration length scale of our 
conformal fluid as a function of intensive fluid parameters. 

In the case of jVA = 4 Yang Mills, it follows from 't Hooft scaling and dimensional analysis that, at 
large A, the effective equilibration length scale is given by Imjp = f{v)/T where v is the dimensionless 
chemical potential conjugate to the conserved charge of the theory. Explicit computation within 
gravity demonstrates that f{v) is of unit order for generic values of v. Consequently, at generic values 
of y, all the considerations of this paper apply only when all fields vary at distances and times that 
are large compared to the local effective temperature. However, as was explained in detail in [42], 
the charged fiuid we study in this paper has an upper bound on v at v^- At this special value of v, 
the black brane becomes extremal. f{v) appears to have a simple zero about nuc- Consequently it 
appears to be possible to scale v to Vc and T simultaneously to zero while keeping Imfp (and all other 
thermodynamic densities) finite. Thus it would naively appear that the long distance field theory 
dynamics should be well described by fluid dynamics in this coordinated extremal limit. It turns out, 
however, that the bulk gravitational solutions described in this paper turn singular in the same limit. 
We are not completely sure how to interpret this fact. It would certainly be interesting to investigate 
this further. 

Note Added : While this draft was in preparation, we became aware of a similar work by Haack 
et.al.(arXiv:0809.2488). 

2. Notations and Conventions 

In this section, we will establish the basic conventions and notations that we will use in the rest of 
the paper. We start with the five-dimensional action^ 



■^We use Latin letters A,B^ {r,v,x,y,z} to denote the bulk indices and /i, € {v,x,y,z} to denote the boundary 
indices. 
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which is a consistent truncation of IIB SUGRA Lagrangian on AdSs x background with a cosmo- 
logical constant A = — 6 and the Chern-Simons parameter k = — 1/(2-\/3) ( [45] - [52]). However, for 
the sake of generality (and to keep track of the effects of the Chern-Simons term), we will work with 
an arbitrary value of k in the following. 

The field equations corresponding to the above action are 



Gab - 6gAB + 2 



FacF^b + -^QabFcdF^^ 



= 

(2.2) 



WBF^'' + Ke^^'^°^FBcFnE=0 

where gAB is the five-dimensional metric, Gab is the five dimensional Einstein tensor. These equations 
admit an AdS-Reisner-Nordstrom black-brane solution 

ds^ = -2u^dx^dr - r'^V{r, m, q) UnU^dx^dx" + r"^ P^^dx^^ dx" 

V3g (2-3) 

where 

r.2 



Uij,dx'^ = -dv; V{r, m,q) = 1 j + -q 



m q 

(2.4) 



with rjni, = diag( — h ++) being the Minkowski- metric. Following the procedure elucidated in [1], we 
shall take this flat black-brane metric as our zeroth order metric/gauge field ansatz and promote the 
parameters u^,m and q to slowly varying fields. 

In the course of our calculations, we will find it convenient to use the following 'rescaled' variables 

''^k' ^^k^' ^' = ^-1 (2-5) 

where R is the radius of the outer horizon, i.e., the largest positive root of the equation V = Q. In 
terms of the rescaled variables, the outer and the inner horizon are given by 

p+ = l and p_ = (g2 + l/4)^^^-l/2 

and the extremality condition p+ = /9_ corresponds to (Q^ = 2,M = 3). We shall assume the black- 
branes and the corresponding fluids to be non-extremal unless otherwise specified - this corresponds 
to the regime < < 2 or < M < 3 which we will assume henceforth. 

Using the flat black-brane solutions with slowly varying velocity,temperature and charge fields, 
our intention is to systematically determine the corrections to the metric and the gauge field in a 
derivative expansion. More precisely, we expand the metric and the gauge field in terms of derivatives 
of velocity, temperature and charge fields of the fluid as 

Qab = 9aI + 9aI +9aI + --- , . 

/l,^_^(0)+ 4(1) , ^(2) , ^ ■ 

(k) (k) 

where and A)^' contain the k-th derivatives of the velocity, temperature and the charge fields 
with 

g'xldx'^dx^ = -2uf,{x)dx^dr - r'^V{r,m{x),q{x)) Ufj,{x)uv{x)dx''dx'" + r'^Pfj,t,{x)dx^dx'' 



...... ^'-'^ 

2r2 



<)rf,M ^ ^u,{x)dx^ 
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In order to solve the Einstein-Maxwell-Chern-Simons system of equations, it is necessary to work 
in a particular gauge for the metric and the gauge fields. Following [1], we choose our gauge to be 



9rr — 



0; gr^ oc ; 



Ar = 0; Tr[(ff(°))-i<?^'=)]=0 



(2.8) 



Further, in order to relate the bulk dynamics to boundary hydrodynamics, it is useful to param- 
eterise the fluid dynamics in the boundary in terms of a 'fluid velocity' u^. In case of relativistic 
fluids with conserved charges, there are two widely used conventions of how the fluid velocity should 
be defined. In this paper, we will work with the Landau frame velocity where the fluid velocity is 
defined with reference to the energy transport. In a more practical sense working in the Landau frame 
amounts to taking the unit time-like eigenvector of the energy-momentum tensor at a point to be the 
fluid velocity at that point. 

Alternatively, one could work in the "Eckart frame" where the fluid velocity is deflned with 
reference to the charge transport where the unit time-like vector along the charged current to be the 
definition of fiuid velocity. Though the later is often the more natural convention in the context of 
charged fluids, we choose to use the Landau's convention for the ease of comparison with the other 
literature. We will leave the conversion to the more natural Eckart frame to the future work. 

In the next section, we will report in some detail the calculations relevant to finding the first order 
solution. 

3. First Order Calculations 

In this section, we present the computation of the metric and the gauge field up to first order in 
derivative expansion, the derivative being taken with respect to the boundary coordinates. We choose 
the boundary coordinates such that = (1, 0, 0, 0) at x/'-. Since our procedure is ultra local therefore 
we intend to solve for the metric and the gauge field at first order about this special point x^. We 
shall then write the result thus obtained in a covariant form which will be valid for arbitrary choice 
of boundary coordinates. 

In order to implement this procedure we require the zeroth order metric and gauge field expanded 
upto first order. For this we recall that the parameters m, q and the velocities {(3i) are functions of the 
boundary coordinates and therefore admit an expansion in terms of the boundary derivatives. These 
parameters expanded up to first order is given by, 



The zeroth order metric expanded about up to first order in boundary coordinates is given by 



m 



q 



mo + aj^a^m^o) + . . . 
qo + x^'d^q^''^ +... 
x^'d^lSi + ... 



(3.1) 



-2x^ d^pf^ dx' dr^2 ^ r^{l - (r)) dx' dv 

(-x>^djM^ 2qoX^dj^\ 2 




(3.2) 
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where mo and qo are related to the mass and charge of the background blackbrane respectively and 
Similarly the zeroth order gauge fields expanded about up to first order is given by, 



2 



\dv- -^x'^Ofj, Pi dx"- 



(3.3) 



Since the background black brane metric preserves an 5*0(3) symmetry ^, the Einstein-Maxwell 
equations separate into equations in scalar, transverse vector and the symmetric traceless transverse 
tensor sectors. This in turn allows us to solve separately for SO{2>) scalar, vector and symmetric 
traceless tensor components of the metric and the gauge field. 

3.1 Scalars Of 50(3) 

The scalar components of first order metric and gauge field perturbations {g'^^^ and A'^^^ respectively) 
are parameterized by the functions h\{r), k\(r) and wiir) as follows ^, 



Y^g^{r)=^T^h,{T), 

i 



(3.4) 



/ii(r) 



4')(r) 



\f^w\ (r) 



Note that g'^\r) and (r) are related to each other by the gauge choice Tr[{g^^^) ^fif^^^] = 0. 



Constraint equations 

We begin by finding the constraint equations that constrain various derivatives velocity,temperature 
and charge that appear in the first order scalar sector. The constraint equations are obtained by taking 
a dot of the Einstein and Maxwell equations with the vector dual to the one form dr. If we denote 
the Einstein and the Maxwell equations by Eab = and Mab = 0, then there are three constraint 
relations. 

Two of them come from Einstein equations. They are given by, 

g^'^E^r + g^^E^y = 0, (3.5) 

and, 



5'''-£^^+5™i;„^ = 0. (3.6) 

^We axe erefering to here to the SO (3) rotational symmetry in the boundary spatial coordinates. 
*here i runs over the boundary spatial coordinates, v is the boundary time coordinate and r is the radial coordinate 
in the bulk 
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and the third constraint relation comes from Maxwell equations and is given by, 



g-'Mr + 9'-''M, = Q . (3.7) 

Equation (3.5) reduces to 

a.mW = ~modif3f . (3.8) 

which is same as the conservation of energy in the boundary at the first order in the derivative 
expansion, i.e., the above equation is identical to the constraint 

d^Tl^o) = • (3-9) 

on the allowed boundary data. 

The second constraint equation (3.6) in scalar sector implies a relation between /ii(r) and fci(r). 

2ai/3f + 12r^hi{r) + 4goWi(r) - mor^h[{r) + Sr'^h[{r) - r^k[{r) - 2qorw[{r) = 0. (3.10) 

The constraint relation coming from Maxwell equation (See Eq. (3.7)) gives 

5^°) = -qo diPf^ . (3.11) 

This equation can be interpreted as the conservation of boundary current density at the first order in 
the derivative expansion. 

Jfo) = 0. (3.12) 

We now proceed to find the scalar part of the metric dual to a fluid configuration which obeys the 
above constraints. 

Dynamical equations and the solutions 

Among the Einstein equations four are 5*0(3) scalars (namely the vv, rv, rr components and the trace 
over the boundary spatial part). Further the r and w-components of the Maxwell equations constitute 
two other equations in this sector. Two specific linear combination of the rr and vv components of 
the Einstein equations constitute the two constraint equations in (3.8). Further, a linear combination 
of the r and t^-components of the Maxwell equations appear as a constraint equation in (3.11). Now 
among the six equations in the scalar sector we can use any three to solve for the unknown functions 
hi{r), ki{r) and Wi{t) and wc must make sure that the solution satisfies the rest. The simplest two 
equations among these dynamical equations are, 

5h[{r)+rh'l{r) =0. (3.13) 

which comes from the rr-component of the Einstein equation and, 

&qoK{r) + w[{r) - rw'({r) = 0. (3.14) 

which comes from the r-componcnts of the Maxwell equation. We intend to use these dynamical 
equations (3.13), (3.14) along with one of the constraint equations in (3.8) to solve for the unknown 
functions hi{r), fci(r) and wi{r). 
Solving (3.13) we get, 

Mr) = ^+Cl, (3.15) 
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where Cj^^ and are constants to be determined. We can set C^^ to zero as it will lead to a non- 
normalizable mode of the metric. We then substitute the solution for /ii(r) from (3.15) into (3.14) 
and solve the resultant equation for wi{r). The solution that we obtain is given by, 

Mr) = Clj' + Cl^-qo^. (3.16) 

Here again C^^ , C^,^ arc constants to be determined. Again C^^ corresponds to a non-normalizable 
mode of the gauge field and therefore can be set to zero. 

Finally plugging in these solutions for hi{r) and Wi{r) into one of the constraint equations in (3.8) 
and then solving the subsequent equation we obtain, 

2 



hir) = Ir^O^pr^ + - 'fcl + r-^-^]cl (3.17) 



Now the constants Cki and C^^ may be absorbed into redefinitions of mass (mo) and charge (qo) 
respectively and hence may be set to zero. Further we can gauge away the constant C^^ by the 
following redefinition of the r coordinate, 

C being a suitably chosen constant. 

Thus we conclude that all the arbitrary constants in this sector can be set to zero and therefore 
our solutions may be summarized as, 

hi{r) = 0, wi{r) = 0, fci(r) = ^r^a,/?^ . (3.18) 
In terms of the first order metric and gauge field this result reduces to, 

i 

gilHr) = ^rd.f3f\ 
9il\r) = 0, 

A(i)(r) = . (3.19) 
Now, we proceed to solving the equations in the vector sector. 
3.2 Vectors Of SO{3) 

The vector components of metric and gauge field g^^^ and A^^^ are parameterized by the functions 
jl^\r) and gl^\r) as follows. 



<'2'W = (?^-?)a!«W 



\ f ^/3(7o\ .(1)/ N , (1)/ X (3.20) 



Now we intend to solve for the functions ij^^r) and gf'^{r). 
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Constraint equations 

The constraint equations in the vector sector comes only from the Einstein equation. So there is only 
one constraint equation in this sector. It is given by 

5'^'-£rt+5™Ki = (3.21) 

which implies 

a^mW = -4moa„/?f \ (3.22) 

These equations also follow from the conservation of boundary stress tensor at first order. We shall 
use this constraint equation to simplify the dynamical equations in the vector sector. 

Dynamical equations and their solutions 

In the vector sector we have two equations from Einstein equations (the ri and ui-components) and 
one from Maxwell equations (the ith-component) ^. 

The dynamical equation obtained from the t;i-component of the Einstein equations is given by, 

- 3my) ^ + 4V3,o.^^MW ^ _ ^ _3,4,^^(o). (3 ^3) 

Also the dynamical equation from the zth-component of the Maxwell equation is given by, 
r [2 {r' - mor' + ql) + (6r^ + 2mor' - 6%^) ^^^^ 



- ^/3gor (r« - mor' + q^) + VSqo {r' - Smor' + 5ql) 

= ^/3(goa„/3f ^ + d^q^'^y - 24ql^rlP, 

where U is defined as 

k = ^ijkdjPk- (3.25) 

Now in order to solve this coupled set of differential equations (3.23) and (3.24) we shall substitute 
gi^\r) obtained from (3.23) into (3.24) and solve the resultant equation for ji^\r). For any function 
jl^\r), using (3.23) g\^\r) may be expressed as. 



9i'\r) = {C,), + ^ I + 4moj W(r) - (^"^^ " ^o) ^ | . (3.26) 



Here (Cg)i is an arbitrary constant. It corresponds to non normalizable mode of the gauge field and 
hence may be set to zero. 



^Note that a linear combination of the ri and ui-components of the Einstein equation appear as the constraint 
equation in (3.22). 
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Substituting this expression for g^^\r) into (3.24) we obtain the foUowing differential equation for 
{35qt + (r^ - 6mo) q'o + Smor^ (3r^ + mo)) "^^''^^"^^ 



r (-ll^o - (5r^ - Wmor^) - niar^ {r'^ + 3mo)) 



dr 



J3„■(l)/„^ 

+ {q^ - mor^) (r^ - mor^ + gg) 



(3.27) 



dr^ 

^ (6\/3go9ig(°V^ + sVsd^pP {5r' - mor' - ql) - 144 r ifql^) 



V3 

The solution to this equation is given by. 



v/3 ;f 6r^go(9,g(°) + 3go9,/3f ^) r mo 



(3.28) 



where again (Cj)i and (C'^)i are arbitrary constants. (C'|)i corresponds to a non-normalizable mode 
of the metric and so is set to zero. (C|)i can be absorbed into a redefinition of the velocities and hence 
is also set to zero. 

.,6 



Here the function F\{-^,^) is given by' 



p I 

where M - 1. 

Substituting this result for j'f'^ {r) into (3.26) we obtain the following expression for gpir), 

= .V3qoRHRlLoir^-I^)) N i'^' + ^^-) ^0 

+<jo(a,:<?(") + 3goa./3f ) - moii^ + mor^) Fi'^'\p, ^^'^^^ 

+6Fi{p,M)Rl (mo-i?cl) r^)] 

where we use the notation f^^'^\a,P) to denote the partial derivative d^^^ f /dd^d(3^ of the function 
/• 

Plugging back jf^{r) and S'-^^(r) back into (3.20) we conclude that the first order metric and 
gauge field in the vector sector is given by, 

q(0) I \/3/f^ggAc „(o) , o„ ^ /l^")^ r ""^ 

ff« - rCvPi — — 4 ^ -wrQa(Oiq' + ^qoOvPi )-ti[ — , -^j 

mor Kg Ho -Kq 

^(i)(^) ^ -y3mogo(9ig^°^ + 3goa./?f Vi''-°^(i^, + IS^^gg/^ 

' 12mQq()r^ 



^Although the expression for ^) is very complicated but it satisfies some identities. One can use those 

rtO -Kg 

identities and simplify the expressions. 
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3.3 Tensors Of S0{3) 

The tensor components of the first order metric is parameterized by the function a\j^ (r) such that, 



(1) 2 (1) 

n..= r ct- ■ . 



(3.32) 



The gauge field does not have any tensor components therefore in this sector there is only one unknown 
function to be determined. 

There are no constraint equations in this sector and the only dynamical equation is obtained from 
the ij-component of the Einstein equation. This equation is given by, 



dr 



c (0) 4 



where cry is given by. 



(3.33) 



(3.34) 



The solution to equation (3.33) obtained by demanding regularity at the future event horizon and 
appropriate normalizability at infinity. The solution is given by. 



.,(1) 



where the function F2 (p, M) is given by, 

F2{p,M), 



p{p'^ +P + 1) 



(p + l)(p4+p2_M + l) 



dp 



(3.35) 



(3.36) 



with M = m/R^ as before. 

Thus the tensor part of the first order metric is determined to be. 



2r 



r mo , 



-5-o-ij-F2(-5-, -^4)- 
no ^0 -n-o 



3.4 The global metric/gauge field at first order 



(3.37) 



In this subsection, we gather the results of our previous sections to write down the entire metric and 
the gauge field accurate upto first order in the derivative expansion. 
We get the metric as 



ds = gAsdx dx 

= -2u^dx''dr - V u^.u^dx^dx" + r'^P^^dx^'dx" 



- 2u^dx'^ r 



2r^ 

dx" + —F2{p, M)af,^dx^'dx'' 



2Uf,dx^ 



mr 



4 L + ^P^VxqF,{p,M) 



(3.38) 



dx" + ... 



dxf" 



where we have defined 
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2 



and 



p4 + p6 j 



(3.40) 



p (p^ + p + 1) 

3.5 The Stress Tensor and Cheirge Current at first order 

In this section, we obtain the stress tensor and the charge current from the metric and the gauge field. 
The stress tensor can be obtained from the extrinsic curvattire after subtraction of the appropriate 
counterterms. We get the first order stress tensor as 

T^v = piVnu + ^u^u^) - 2j7C7^^ + . . . (3.41) 

where the fluid pressure p and the viscosity ry arc given by the expressions 

R"^ s 
p = — — ; r? = — — = — (3.42) 

where s is the entropy density of the fluid obtained from the Bekenstein formula. 
To obtain the charge current, we use 

= lim =nu^-D P;V,n + ^1^ + ... (3.43) 

where the charge density n, the diffusion constant 5) and an additional transport coefBcient ^ for the 
fluid under consideration are given by 



IGTrGg AMR ' ^ IGTrGsm 

We note that when the bulk Chern-Simons coupling k is non-zero, apart from the conventional 
diffusive transport, there is an additional non-dissipative contribution to the charge current which is 
proportional to the vorticity of the fluid. To the extent we know of, this is a hitherto unknown effect 
in the hydrodynamics which is exhibited by the conformal fluid made of A/" = 4 SYM matter. It would 
be interesting to find a direct boundary reasoning that would lead to the presence of such a term - 
however, as of yet, we do not have such an explanation and we hope to return to this issue in future. 

The presence of such an effect was indirectly observed by the authors of [42] where they noted a 
discrepancy between the thermodynamics of charged rotating AdS blackholes and the fluid dynamical 
prediction with the third term in the charge current absent. We have verifled that this discrepancy 
is resolved once we take into account the effect of the third term in the thermodynamics of the 
rotating J\f = A SYM fluid. In fact, one could go further and compare the first order metric that we 
have obtained with rotating blackhole metrics written in an appropriate gauge. We have done this 
comparison up to first order and we find that the metrics agree up to that order. 
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4. Second Order Hydrodynamics 



In this section we will find out the metric, stress tensor and charge current at second order in derivative 
expansion. We will follow the same procedure as in [1] but in presence charge parameter q. 
We will consider the following metric and gauge field perturbations in second order, 

g'^^^^dx^dx^ = -3h2{r)dvdr + r'^h2{r)dx'dxi + '^^dv'^ + Ur'^j^dydx' + r^a^fdx'dx^ (4.1) 

and 

^(2) = fr^gPir)dx^ . (4.2) 

Here we have used a little different parameterizations (from first order) for metric and gauge field 
perturbations in the vector sector. We found that the corresponding dynamical equations for jf^\r) 
and gl'^\r) can be written in more tractable form (wo will see later). 

Like neutral black brane case, here also we will list all the source terms (second order in derivative 
expansion) which will appear on the right hand side of the constraint dynamical equations in scalar, 
vector and tensor sectors. These source terms are built out of second derivatives of m, q and /3 or 
square of first derivatives of these three fields. We can group these source terms according to their 
transformation properties under SO(3) group. A complete list has been provided in table 4. li and 
aij are given by, 

h = £ijkdj(3k , aij = ^{dil3j + dj(3k) - ^Sijdk^k ■ (4.3) 

In this table we have already employed the first order conservation relations i.e. equation 3.9 and 
3.10. Using these two relations we have eliminated the first derivatives of m and q. At second order 
in derivative expansion we have to use the relations, 

d^d^T^,"^ = (4.4) 

and 

d^d^J^„^ = . (4.5) 

These two equations imply some relations between the second order source terms listed in the 
table. 



SI 

S2 



^-^STl 



16™ST3-^ST4+ll^ST5 



9 



--r-S3 

4m 



-STA 
2 



4ST1 

QSl = q (-ST1 - S2 + ST3) - QS5 
4C 
"9 



3 

ST5 



VI 40,,, 4,,^ 56,,^, 4,,_ 8,,^„ 

— = V4 - -V5 H VTl + -VT2 + -VT3 

m 9 9 3 3 3 

in 1 1 ^ 

V2 = — V4 + -V5 - -VTl + -VT2 - -VT3 
9 9 3 6 3 



(4.6) 
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1 of SO{3) 


3 of 50(3) 


5 of 50(3) 


01 

oi = 




Vlj = didym 


rp-| 


OiOjm — 3 


00 

bZ = 


fi fi ^ 

OvOiPi 




rpn 

— 




bo = 


am 


V3i = dyli 


rpq 

LOij — 




oil 





V4i = pjaji - d^di 


rprpi 


= OyPiOyPj - 3 bil djj 


Qrnrt 
biz 


7 Q /a_ 


V5i = d'^pi 


rprpo 

TT2ij 


— l^iOvPj) — 3 bii Ojj 


bio 




VTli = ^{dypi){djp^) 


' 1 " 1 10 
i i 6ij 


= ^^ki{i<JvP Oj)P +3 oiz 


b ± 4 


— H I 


VT2i = -djk P dyl3^ 


i i4ij 


— OkP <7ij 


1 


— cFij a 


VT3i = Uij dyP^ 


± i Oij 


— H 3 b±4 Ojj 




— rfin 


VT4i = h djpi 


1 i Kiij 


_ _.. _fe 1 CTc; A. . 
— CTife (Tj — 3 010 Oij 




= didiQ 


VT5i = Gij P 


TT7. . 
i i 




QS3 


= kdiQ 


QVli = diqdkP'' 




= didjq - i QS2 5ij 


QS4 


= {diqf 


QV2, = diqdkf3'' 




= d{iqlj) — ^ QS3 5ij 


QS5 




QV3j = Cijkdjlk 


QT3,^. 


= d{iqdj)q- \ QS4 5ij 






QV4i = aijdjq 


QT4,^- 


= dij^qdyfij) - \ QS5 (5jj 






QV5j = eijkdy(3jdkq 




= £{ikm9kq O-mj) 



Table 1: An exhaustive list of two derivative terms in made up from the mass, charge and velocity fields. In 
order to present the results economically, wc have dropped the superscript on the velocities /3i charge q and the 
mass m, leaving it implicit that these expressions are only valid at second order in the derivative expansion. 



V3 = — VT4 + VT5 

3 

QVf = q (^y V4 + ^(VT2 + 2VT1 + 2VT3) + iv5 
+QV2 + i (^2QV4 + QV3 + |qV2 
Tl = -4m ( T3 + ^TTS - 4TT1 + ^TT4 + TT6 



(4.7) 
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With these relation between the source terms we will now solve the Einstein equations and Maxwell 
equations to find out the constraint and dynamical equations at second order in derivative expansion. 
We will do the analysis sector by sector. 

4.1 Scalar Sector 

Our parametrization for metric and gauge field are following, 

J29i^\r)=3r'h,{r), 



Scalar Constraint Equations 



e){r)=-\h2{r) 

V^W2{r) 



Mr) 

(4.8) 



2r2 



As we have already explained, there are three constraint equations. First two come from Einstein 
equations (Eq. 3.5 and 3.5) and the third one comes from Maxwell equations (Eq. 3.7). The first 
constrain from Einstein equations gives, 

a„m(i) = \rI ST5 (4.9) 
Second constraint implies relation between k2{r) and /i2(r). This constraint equation is given by, 

k'^ir) = Sk (4.10) 

where Sk is given in appendix A. 

The constraint relation coming from Maxwell equations is given by, 

_ 3go(flo^ + mo) S3 (V + mo) 6V3go^^ 
lomo^iio 4mo-Ko tio 

+ ^/^K-11^o^)qS5 _ ^QS3 ^QS4 

SmoRo mo 4mo-Ro 

9qo iiRf)' + mo) , ^ 

_yu_V — u u^g^^ 

4moi?o ^ ' 

Scalar Dynamical Equations 

Scalar Dynamical Equation is given by {Err = 0), 

f>{h\{r)+r{h'\{r) = SH. (4.12) 

The source term Sh is given in appendix A. 

The second dynamical equation in this sector, which is M(r) = 0, is given by^, 

dw2{r) (Pw2(r) dh2{r) 2 1 



''There is a relative negative sign in the r.h.s. with respect to first order equation and this is because of different 
parametrization. 
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The source term Suir) is given in appendix A. 

The source terms have the same large r behavior as uncharged case because the charge dependent 
terms (leading) arc more suppressed than that of charge independent terms. So one can follow the 
same procedure to find the solution for h2{r) and fc2(r). 

1 f°° fjry f°° 1 

^2(0 = -4^5,- + 1 dyy\S,iy) - -^Sn (4.14) 

and 

fe(r) = - dxiSkix) - xSn (4.15) 

4.2 Vector Sector 

Vector Sector Constraint Relation 



The constraint relation coming from Einstein equations (3.21) is give by, 

+ ^S!A^QV4. - <^'f -'■='"-f) vT3. . (4.16) 
mo — 3i?o 3 (mo — 3i?o j 



Vector Sector Dynamical Equations 



There are two vector dynamical equations. The first equation which comes from Einstein equation 
is given by, 

3r [5rqoM'\r) + r^9o^ + 5^ + r^^^j = Si^\r) (4.17) 

where Si^'' is the source; terms given in the appendix B. The second dynamical equation comes from 
Maxwell equation and is given by, 

5V3 (7r^ - 3mor2 + go') Q^ir) 

12V39o^V^ + V3r (I3r« - 9mor^ + 7qo') 

dr ^ ' dr 

x/3r^ (r*^ - mor^ + q^) ^^^^ = -2S^^\r) (4.18) 

where Sy is other source terms given in the appendix B. Prom this above two equations one can 
replace g\ ^ (r) and write a differential equation for (^) • 

3V3 {-2.hr^ + 15mor2 + W) - 3V3r (l3r« - Qmor' + 7go') ^^4^ 

dr dr^ 

-3x/3r2 (,6_mor2 + <^o2)^!,g(!l 
= - f2V3r (mo - 3r^) S^^\r) + H^rSf\r) -^J>{r^- m^r" + go') ^^^) (4-19) 
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This equation can be written as 
d f Id 



dr \ SmoRo dr 



= 'S'source (4.20) 

4.3 Charge Current 

To find the charge current density we need to solve the coupled dynamical equations (Eq.4.17 and 
Eq.4.18) in this sector and wc have to consider the asymptotic limit. The charge current at second 
order in derivative expansion is given by, 

j(2) = LJh^ . (4.21) 

Now all the computation can be done setting Rq = 1. In the final solution the correct power of Rq 
can be restored using Weyl covariance of the current. The solution to the Eq. 4.20 can be written as 
follows, 

where M, Q and p are given by, M=^ , Q = P — ^ ■ "^^^ function /i(x) is given by, 

,.5 



hix) = _ ^ ^ 2 (^Ms(i)x4 - 8Ms(2)a; - 4Ms(3) log(a;) + SM J y H{y) dy + 



(a;6 - Ma;2 + Q'^f 
where H{y) is given by, 



(4.23) 



S^oru.{x)dx (4.24) 

and the constant C is given by, 

C = -S(i) (2M - 3) + 2s(2) (5M - 3) + S(3) (M - 3) + B (4.25) 

where, 

/oo roc 

SNorm{x)dx + 4M J x'^ S^orn,{x)dx . (4.26) 

We choose this constant C in such a way to cancel the divergence (in this case logarithmic) at the 
horizon Rq = 1 for the function jl'^\r). The function ^Norm is the normalized part of the source terms 
where we have subtracted out the asymptotic divergent parts. 

Snotui = SsouTce ~ r ^jg- . (4.27) 

The solution to the gauge field perturbation is given by, 
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where, 

SilM = Si'^ - '-^ . (4.29) 

The charge current can now be conveniently wriiten down in a manifestly Weyl-covariant form. 
The charge current is given by, 

^f^ = -^EC^^» (4-30) 



where, 



Li — —Ho 




4M 

2VM^ (3(M - 1)(2M - 3)k2 - M^) 

C3 = 

C4 = 
C5 = 

Ce = i?o "° (4.31) 

Here -B, is the coefficient of Wj in B (defined in (4.26)) and the Weyl-covariant vector combinations 
W, are given by 

5V4 5V5 5VT1 5VT2 11VT3 , 

^^ = ^+— +^ — 12 — ^=^;^^-^^ 

5V4 V5 VT VT3 x 

H^. = — -^-VTi-- + ^ = p;W 

OV3 3VT2 
1^3 = ^ + QV4 + + 3qVT3 = -P^u^V^Vxq 

W4 = QV4 + 3qVT3 = a^^Vxq 
Wg = QV3 + 3qVT2 = w^^Pa? 

We = VT5 = u^Vxl^ (4.32) 
where we have written the vectors in terms of the Weyl-covariant derivative introduced in [6]. 
4.4 Tensors Of 5*0(3) 

We now consider the tensor modes at second order. Following the first order calculations we parametrize 
the traceless symmetric tensor components of the second order metric by the function Tij{r) such that, 

i?l?=r-%- (4.33) 

The ij-component of the Einstein equation gives the dynamical equation for which is similar 
to (3.33). However the source term of the differential equation is modified in the second order. This 
equation is given by, 
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(r^ - mor^ + ql) r'l^ (r)r - (-5r*^ + mor^ + go') r/,- (r") = T^,- (r). (4.34) 

where, 



Ti, (r) = Ci(r) (^TTly + irr4i,- + T3y^ + C2(r) TTh^ + C3(r) TTe^ + C4(r) TTTi,- 

+ C5(r) (QT3y + 6goQT4i,- + Qg^TTly) + C6(r) (QTl^,- + 12goTTly + 8QT4i,) 
+ C7(r) (QT2y +3goTT2y) + C8(r) (4TT2y + 2T2y - TT3y) . 

(4.35) 

with the coefficient functions, 



Ci(r) 


_ (' 




3r , 






C2(r) 







(r + i?o)(r4 + i?2r2 + i?4_^^) 



144K2g6 36^2^4 36^2^4 27rg2Fi (-^, ^) 



(1,0) I r mp 
Rq ' -Rq 



mor^° mgr^ mor^ i?5 

q^2p(l,0) / j_ rna\ 9 

-^1 l^flo'flfjgp _ 2g| _ mo _ 

= - Rtir + Rp)ir^ + Rlr^+Rt-mo) ' ^^^^^ + " + + ^^^^ 

-3r(r + i?o) (r^ + i?gr^ + i?^ - -o) ^) i?^ + 9 (i?o^ - mo) F, ^) iio 

.3.(.^-mo).r(£,^))) 



C4(r) 



2Rl 



-»('^-'^))H'"fei) (i-i))) 

- (" - ^- (i' i) - ^^•"» - (i' i) 
(w (^, ^) «j + (^, ^)) ^ 

.2 d8 _ . o4 _L p(l,0) / mo \ ' _^ pio _ 1 n^. U6\ 17(1.0) / ™0 ' 



+3r2 (3ii« - 4moi?^ + ml) [jT^^^J + (8^o° - WruoRl) F^ (^^, ^ 



(4.36) 
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Bl 



+2r3 (-r^ + i?6 + mo(r - i?o)(r + i?o)) 



Cs{r) = 



(2,0) j r uiQ 

^—1 -r -11-0 -r "'01.' ~ -"-OA' "r -"-o;; 

2mor^ 

(4.37) 



The solution to (4.34) which is regular at the outer horizon and normalizable at infinity is given 

by, 

-« = -f <«(^5.TieTF)£''C2CT«((). (4.38) 

For the computation of the second order stress tensor of the boundary fluid we require the large r 
behavior of only. Using (4.38) and the coefficient functions in (4.36) we obtain such an asymptotic 
expansion for which is given by, 



+ ^ {^TT^ ™ij + ^TT7 TTlij + Nttq TTGij + AArs (rTlij + ^TT4ij + T3ij^ (4.39) 

+J^QT2 {QT2ij + 3qoTT2ij) + Mqti {QThj + UqoTTUj + 8QT4ij) 

+Mtts (4TT2ij- + 2T2ij - TT3ij + Mqts {QT3ij + QqoQTAij + Qg^TTli^))) . 



where, 



(mo - Rt) {2lRt - 672 (mp - Rf) + 53mo) 

" 224moi?;^ 

A/-,,, =1 (^,(1, ^) _ [p^Fr\p, ^) + 



1 / / \ r 2(l + p + p2 + 2p4)-^(l + p)^ 

^--75 H *(.^^.(..^)-.)+/ 

-"0/ -"0 ^-"0 



(4.40) 



i?0 

^^^^^^ (^A 



r,/,^, mo, ^ ^oo p(l + p + p^)(2 + 3p-^)-^(l + p) - 
A (l+p)(,^+p^ + l-^) 



+ 



/ 3mo 3 3i?o 
V56i?o ~ 7^ 8^ 
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-ftf, -Kn 



TT3 



(4.41) 



4moi?g 

1 r°° 

J^QTs. =77^ drr C5(r) 



The leading term of this asymptotic expansion gives divergent contributions to the stress tensor 
which are cancelled by divergence arising from the expansion of 5^°^ + g^^^ up to second order. The 
subleading term in (4.39) yield the second order stress tensor. 

4.5 Stress Tensor 

The stress tensor is obtained by the following formula, 



= -2 lim (Kj^ - Sj^) , (4.42) 

r — *oo 



where Kj^ is the extrinsic curvature normal to the constant r surface. If we evaluate this expressioin 
on our solution the second order stress tensor is obtained to be proportional to the coefficient of the 
^ term in (4.39) multiplied by a factor of 4. 

The stress tensor so obtained can be shown to be Weyl-covariant using the following relations 

u^T>x(Ti,^ = TTl + \tT4: + T3. 

cT^xcTx. + Ip^-'ct'^'^ct^p = TT6. 
-2 {w'^a^-' + LO^^a^^") =TT7 

I [P^PI' + P^P^ - Ip"^'P^.^ V^qVpq = QP3 + 6ggT4 + Og^TTl (4.43) 
\ [P^P^ + P^P^ - ^P^^'P^.) V^V^q - 3 q u^V^a^, 

-3 q (^ui'xux^ + ip^''w"'^w„/3^ = QTl + SQTA + UqTTl 

'D^l^ + VJ^ = ATT2 + 2T2 - TT3. 

\ {liJ>uq + lu'D^q) = QT2 + 3q TT2. 

These relations can be used to express the energy-momentum tensor in the manifestly Weyl-covariant 
way. 

5. Discussion 

In this paper, we have computed the metric dual to a fluid with a globally conserved charge and used 
that to find the energy-momentum tensor and the charge current in arbitrary fluid configurations to 
second order in the boundary derivative expansion. 
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Note that the corresponding construction of the bulk metric dual to an uncharged fluid flow was 
characterized by a great deal of universality. This universality had its origin in the fact that every two 

derivative theory of gravity (interacting with other fields) that admits AdSs as a solution, also admits a 
consistent truncation to the equations of Einstein gravity with a negative cosmological constant. This 
universality does not extend to Einstein Maxwell system. The possibility of extending the Einstein 
Maxwell system by a Chern Simon's term of arbitrary coeflicient, accounted for in this paper, is an 
illustration of the reduced universality of our calculations. 

We have seen that a nonzero value for the coefficient of the Chern-Simons in the bulk leads to an 
interesting dual hydrodynamic effect (recall that this coefflcient is indeed nonzero in strongly coupled 
M = A Yang Mills ). At first order in the derivative expansion we find that the charge current 
has a term proportional to P = e^^'^^'^u^jV yU\ in addition to the more familiar Fick type diffusive 
term. It would be interesting to reproduce this term from a computation directly in the boundary. 
In more general terms, it would be interesting to gain better intuition for effect induced by the bulk 
Chern-Simons term on boundary dynamics. 

We have already remarked on the fact that the pseudo-tensor terms appearing in the stress tensor 
and the charge current solve an old puzzle raised in [42] regarding the fluid-gravity correspondence 
in large rotating AdS charged black holes. The discrepancy found by the authors of [42] between the 
fluid dynamical predictions and the known thermodynamics of rotating blackholes is basically resolved, 
once the qualitatively new effects due to the bulk Chern-Simons interaction on hydrodynamics is taken 
into account. We will reserve a deeper analysis of this issue along with a detailed comparison of our 
second order fluid dynamical metric and gauge fleld with charged black hole solutions to future work. 

Though we have not yet worked out explicitly the position of the event horizon in our metric 

solutions, it is plausible that the analysis of [2] can be easily extended to the case of metrics with 
a global conserved charge(at least for the non-extremal case). This expectation, however has to be 
confirmed by an explicit computation. Further, it would be interesting to derive an entropy current 
for the charged fluid following the proposal outlined in [2]. 

In the bulk of our work we have refrained from commenting on the fluids near extremality. The 
preliminary analysis in [42] suggests that hydrodynamics would be a valid description at least for some 
class of extremal solutions. Unfortunately, we have not been able to shed more light on this issue in 
the present work - the perturbative metric and the gauge field we have obtained are ill-behaved at 
the horizon in the near-extremal limit. However, it is not clear what this behavior signifies. This is 
indeed one of the most pressing questions opened up by the present work and it would be interesting 
to understand the extremal limit of our solutions more clearly. 
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A. Source Terms in Scalar Sector: Second Order 

There are three source terms in scalar sector Sk{r), Sh{r) and SM{r). They are quite complicated 
functions. Here we tried to put the explicit form of these source terms as much as possible. 

r (-2flo« + 54goVFl ^) R, + 9goV^Fl(^-°) ^ 



where, 



6moi?o 

- -!—{Ih + Ai+Bi+Ci+Di+Ei) ST5 



S3 



zr Zr / ~ ~ \ 

- — ST3+ — (^i(r)STl + ^2(r)QST5j 

2r3 3qO(6i?oFl(^,i|)+rFl(^'°) 

3V3q03 {mp - r^) . (si^pFlt^-") ^) + rFl^^-") (^, ^g)) 

Ampr^R^ 

V^qO^ (mo - r^) n (siJoFl^^^") f ) + rFl(^.°) (i,, f )) 



ST2 



At = 
Bi = 
Ci = 

Di = 



Ampr'-^RQ 

qp'Rp{r^ + Rpr + Rp''f 
(r + Rpf (r4 + Rp\^ + Rp^ - mpf 
mpr^Rg {r^ +Ror + Ro^f 
(r + i?o)2 (r4 + i?oV2 + i?o* - mo)' 
r^J^o (r^+i?or + flo^)' 
(r + i?o)2 (r4 + Rp\^ + Rp^ - mpf 
Amp {r^ + R,r + Rp')F2U,^ 



QST3 (A.l) 



(r + Rp) (r4 + /?oV2 + Rp^ - m.p) 
12r^r^ + Rpr + Rp^)F2{^,f^) 



(r + i?o) (r4 + i?o V2 + i?, 



-mp 



and the functions Ai(r) are complicated functions. Hence we are not giving the expressions here. But 
these functions have the correct asymptotic behavior. For example for large r, Ai{r) — ^• 
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The second source term is given by, 



where, 



and 



Sh = ^{A + B)ST5 



4^2 



^4 _ 36qn k 

+ . ST4 



i?0 



18 



+ ^ ^ " "q^ ^ " " ST2 

H —g (^b6 

morRo 

(5i^oFl(^-°) +rFl(--°) (i^^l^))' ^,^ 



p 18 
-n-O 



-QS5 (A.3) 



(r + i?o)' {r^ + Ro^r^ + Ro" - mo) 



mo 



(A.4) 



2r (mo (4r3 + Siior^ + &Ro^r + SR^^) - 3Ro^ (r^ + Rqt + i?o^)^) F2 (-^, ^^^^ 
Ro{r + i?o)2 (r4 + Ro\^ + Rq" - mo) 
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'he source term Sm {r) is given by, 



iV^^Rl (mo - Rt) - 12r (mo - R^) Fl (i^, ^)) 

8mor4i?5 ' 

V3(i?§-12r(mo-i?^)Fl(^,^)) 



-QST2 



4V3Vi?^ (mo - i?^) (-r« + i?g + mo (r^ - F2 ^) F2(i.°) f) 
+ 

+ V3{Rl- r^) ^Rlimo-Rl) (r^ + R^ + R^ - mo) 

(36K2i?i2 _ 24too (r^ + 3i?g) K^ije + ^2 / 4 _^ 24i?gKV2 + 36R^k^)) — ^ST4 

2mor^* 



+ 



9V3^R^{mo-Rt){{mo - SR^) Rl' + r (mo - i?^) (24mo (mo - i?^) FlfO'^) (^-^, ^ j i?^ 
+ 3F1 (^-^, ((I3i?^ + 5mo) R^ + 2r^ {3R^ - mo) ((-21i?^ + 21moi?^ - lOmor^iJo) 

Fl^^'°^ ^) + r {-3Rt + 3moi?^ - 2mor^) Fl(^^°) ^) ))i?o 

+ r(i?oFl(i'0) (^-^, ((Umo - 21/?^) + 15r2(r - i?o)(r + i?o) (r* + R^ + Rt - mo) 
(3i?^ - mo) Fl(i'") (^-^, I?) ) - ^ (3i?o - mo) {R'o + 3r^ {-r^ + i?^ + mo(r - i?o)(r + Ro)) 
Fl<^^°^f^,^))Fl^-^f^,^)))) 



Ro'RtJ' \Ro' RtJ"'r^RhHSRt-mo) 

- 9 (mo - -Ro) ^((-5/?^ + 5moi?o - 4mor2) Rq + 2r^(12Rl'^ - 24moi?o + (9r* + Umo) 

Rl + UmlR^ - mor^ (9r^ + Umo) RI + lOmor^ (mo - r^))Yl^^^''^ Ro 

\Ho Ho J 

+ 2r' ((3r^ + mo) R^ - mo (3r* + mo) RI + 2mor^ (mo - r')) Fl^^.o) (^,^^' 



Ro' R^J'2mor^Rl 
V3q0(24mo (mo - R^) Fl'"'^) , ^) R^ + 3F1 , ^) ((E^ + 9mo) 

R^ - 2r' (63< + ml {21RI - lOr') + 6mo (Sr^i?^ - Mi?^)) Fl^^'") ^) Ro 



i?0 ' Rq 

+ 2r' (-9< - 6mo (r^ - 2Eg) i?4 ^ ^2 (2,2 _ 3^2^) ^1(2.0) (^^^ 
+ r(r (mo - 3ii^) Fl(^'°) ^) R^ + 15r\3 (r« - R^) R^ + ml (r^ - R^) 



Ro R'o 



- mo (r« + 3Ry - 4i?^))Fl(i'°) ^) R^ + Fl(i'°) ^) {Rl (llmo - 21i?^) 

- 3r3 (3 (i?^ - r<5) i?^ + (i?^ - r^) + mo (r<5 + 3Ry - AR^)) 

- 3 (mo - i?;|) K((-5i?^ + 5moi?o - 4mor2) R^ + 2r^{3 {3r^ + AR^) R^ 

+ ml (lOr^ - URy + UR^) - mo (lOr« + dR^ - UR^r' + 24i?«))Fl(i'°) (^-^, Ro 

+ 2.^ (3.^..^ H- m^ (2.^ - R^) + mo (-2.^^^i.^.^ + R^)) fK-) )^-^ 

+ v/3(5i?i5 (mo - 3E^) - 12r (mo - i?^) (-24mo (mo - R^) Fl^°''^ (^-^, i?g 

+ r.Fl(i'°) ( ^. f217?i2 _ iimni?^ - ISr^fr - i?n)(r + i?n) Cr^ + iJ^r^ + Ri - mo) 



where, 

func{r) = {3Rt - mo) Fl^'''^ ^) )i?o + 2F1 ^) (Sr^ (SR^ - mo) ((21i?o^ - 21moiio' 

+ mmor'Ro)Fl^''^ {w,^^) + ' ^^^^ " 3moi?§ + 2mor^) Fl^^'") ^) ) 

- R^ {9Rt + llmo))i?o + (3i?^ - mo) {R^ + 3r^ (-r*^ + i?o + mo{r - i?o)(r + i?o)) 

B. Source Terms in Vector Sector: Second Order 

The Source Terms 



" ~ 9r2(r + iio) (r* + i?oV2 + i?o^ - mo) 

(2r5 + 2iior-^ + 2Ro'^r^ - SRo^r'^ - SRo^r - 3i?o^ + 3mo(r + Rq)) V5X 

9r2(r + Rq) {r^ + Ra'^r^ + iio* - mo) 
(2r5 + 2iior-^ + 2Ro^r^ - 3i?o - 3i?o V + 3mor - 3i?o'^ + 3moi?o) VTIX 

3r2(r + Ro) {r^ + Rq^^ + Rq^ - mo) 

50^ (Sr (r2 + Ror + R^^f + mo(r + 2i?o)) VT5Xk 

4\/3mor2(r + Rq)"^ (r* + i?oV2 + Rq- - mof 

TIA + TIB + T2,,^^ 

+ VT3 

TD 

TNIA + TNIB + TiV2 ^, , 
+ QV4 

5r (r2 + i?or + i?o^) 1 ^ 

(B.l) 



6(r + i?o) (r4 + i?o V2 + R^" - mo) 2r2 
3go^QV3X (6F1 f^) R^' + r (9i?oFl(i.°) + rFl(^'°) |^))) 



where. 



Ro^ 



2rRo 



3 (r^ + Ror + i?o^) (Sr*^ + 23i?oV^ - ^Ro^) Ro^ 

3mo^(r + Rof + mo^ (Sr^ + WRor^ - 26i?oV* - 19i?oV^ - 21i?oV2 - 30i?oV - 15i?o^) 
mo {8r^° + 16Ror^ + 24Ro^r^ + 75i?oV^ + 141i?o*r^ + 207i?oV^) 

mo (l09i?oV^ + 32i?oV^ - 45i?oV^ - 42i?oV - 21i?o^°) (B.2) 
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TIB 

-WSqo^r^Ro = ^ + + ^"')' + ^iio') i?o^ + mo^ (Tr^ + 14i?or2 + 12i?oV + 6i?o^) 
- 3mo (r^ + 2i?or^ + 3i?oV^ + 4i?oV'^ + 12i?o^r3 + 20Ro^r'^ + 16Ro^r + 8Ro^) 



T2 = -lOSgoV^ (r^ + 2i?or2 + 2iioV + i?o^) 

(3 (r^ + Ih'r' + Ih') Ro' + mo' - mo (r^ + Ih'r' + 4^)) Fl^^'") (^,— 

(B.4) 

TD = 6r^Ro^{r + Rof (r^ + i^o'r^ + i?o^ - mof (3i?o* - mo) (B.5) 

= i?o'^ + 2r (9 (r^ + iior + Ro'f (r^ + 2iio^) i?o^ + mo^ (7r^ + 14i?or' + 12i?oV + 6i?o^)) 
- 2r (3mo (r^ + 2i?or*^ + 3iioV^ + 4iio^r-^ + l2Ro\^ + 20i?oV^ + 16i?oV + 8i?o^)) 

i?o 



Fl(^.^)^o (B.6) 



TA/'IB = &qor^ (r^ + 2i?oT-^ + 2i?oV + i?o^) 

(3 (r^ + Ro\' + Ro') Ro' + mo' - mo (r^ + Ro'r' + 4Ro')) Fl^^'") ^] 

\Ro Ro J 

(B.7) 

TN2 = qorRo^ (mo(2r + Rq) + 3rR^ {r^ + 2Ror' + 3i?oV + 2i?o^)) (B.8) 

TD = Ro^r + Rof {mo - 3iio^) {r" + Ro'r' + Ro^ - mo)"" (B.9) 
The second source term is given by, 



^ ^ jRo' (mo - Ro^) {24k'' Ro^ + mo (r^ - 24i?o'K2^ ) 

VSmor^ 

VsJRo' {mo - Ro^) {24k'' Ro^ + mo {r' - 24i?o^K')) 

+ ^ ^ — s ^ ^VTl 

mor° 

^ mor6(r + Ro) (r^ + Ro'r^ + Ro^ - mo) 
QV3 



2mor^Ro^ 



+ ^^VT2 

2mor6i?o 

.VT3 



2mor^Ro^ 

S5 



rRo\r + Ro) {r^ + Ro'r' + Ro^ - mo) 



QV4 



+ -F- 7 o -A rV4 (B.IO) 

^/3{r + Ro) (r* + R^'r' + Ro^ - mo) 
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where, 

SI = 3i?o (r^Ro {r {^r^ + '^R^r^ + iRo'^r^ + Ra^r"^ + Ro'^r + i?o^) - mo (Sr^ + ZRqt + 2Ro^)) 

- 8TOo(r + iio)(r^ + i?oV + i?o^-mo))511F2f^,^') (B.ll) 

V-n-o no J 



S2 = V3{{2Ak^R^^ - 24moK'i?o' + mor^) Rq^ + Smor^ {Rq^ - mo) F1 — , -i^ i?o 

V-Kq i?o 



+6mor4 (mo - Ro^) Fl(^'°) ( 



r mO 



(B.12) 



5*3 = V3^JRo^ (mo - i^o^jlS (24K2i?o*^' - 24moK2i?o^ + mor^) i?o^ 
- aemor^ (mo - Ro^) Fl (^-^, i?o 

+ 36mor4(mo-i?o')Fl(i-") (;^,;^)) (B.13) 



S4 = V3^i?oMmo-i?o')(6moFl(2.°) (^-^, r" + 30moi?o"Fl(i>°) 

+ (6mo (r^ + RQr + i?o^) (-i?o^ + 12r (mo - iio*) Fl f ^, ^) Rq^ 

\Ro Rq J 



^ V^o^RoV ^ (r + i?o) (r^ + iio V + i?o^ - mo) 

+iio® (24K2i?o^ + mo (r^ - 24i?o^K^))) (B.14) 

55 = V3{r{12Ro (r^ + i?or + i?o') (mo - Ro^) Fl f 

\Ro Ro J 

+ r{r + Ro) (r^ + Ro'r^ + Ro' - mo) (5i?oFl(^'") (^-^, 

+ ,F1(2'") f ^, ^) )) - iJo*^ (r^ + iJor + i?o^)) (B.15) 

\Ro Ro J 

56 = 5((6 (^^V (mo--Ro^) - \l moRo" - r/^ (12i?o (r' + -Ror + -Ro^) (mo - i?o^) Fl (^-^, 
+ r(r + i?o)(r^ + i?oV + i?o^-mo) (siJoFl^^-^) (i'^) +^^1^''°^ 

- (24(r + i?o) (r^ - i?or + i?o^) (r^ + i?or + i?o^) ^iio^ (mo - V)/«^iio'' 

+ mo2(r + Ro)^Ro^ (mo - i?o^) (24i?o'K' - r^) 

+ mo{^jRo^ (mo-i?o^)(r2 (r^ + i?or + i?o^) (Tr^ + i?o^) 

- 24i?o'(r + Ro) {r" + Ro^r^ + 2Ro^) k^) - 6r^ (r^ + i?or + iio^) moRo" - Ro^))/mor'') (B.16) 
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